Measuring parameters of massive black hole binaries with partially aligned spins 
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The future space-based gravitational wave detector LISA will be able to measure parameters 
of coalescing massive black hole binaries, often to extremely high accuracy. Previous work has 
demonstrated that the black hole spins can have a strong impact on the accuracy of parameter 
measurement. Relativistic spin-induced precession modulates the waveform in a manner which can 
break degeneracies between parameters, in principle significantly improving how well they are mea- 
sured. Recent studies have indicated, however, that spin precession may be weak for an important 
subset of astrophysical binary black holes: those in which the spins are aligned due to interactions 
with gas. In this paper, we examine how well a binary's parameters can be measured when its 
spins are partially aligned and compare results using waveforms that include higher post-Newtonian 
harmonics to those that are truncated at leading quadrupole order. We find that the weakened 
precession can substantially degrade parameter estimation. This degradation is particularly dev- 
astating for the extrinsic parameters sky position and distance. Absent higher harmonics, LISA 
typically localizes the sky position of a nearly aligned binary a factor of ~ 6 less accurately than for 
one in which the spin orientations are random. Our knowledge of a source's sky position will thus 
be worst for the gas-rich systems which are most likely to produce electromagnetic counterparts. 
Fortunately, higher harmonics of the waveform can make up for this degradation. By including har- 
monics beyond the quadrupole in our waveform model, we find that the accuracy with which most of 
the binary's parameters are measured can be substantially improved. In some cases, parameters can 
be measured as well in partially aligned binaries as they can be when the binary spins are random. 
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I. INTRODUCTION 

The coalescence of massive black hole binaries is a pri- 
mary source for the future space-based gravitational wave 
(GW) detector LISA.^ LISA will be able to detect such 
sources with extremely high signal-to-noise ratio (SNR) 
> 1000 at low redshift {z ^ 1), as well as with moderate, 
but still reasonable, SNR 10 at extremely high redshift 
{z ~ 20) [T] . Estimated event rates for these sources vary 
widely based on formation scenarios but tend to predict 
roughly tens of sources per year, with ~ 1 as a pessimistic 
estimate and ^ 100 as an optimistic one [5]. (The actual 
detection rate will, of course, tell us much about the for- 
mation of black hole binaries and the growth of massive 
black holes in the universe.) 

While the detection of gravitational waves from these 
sources will certainly be interesting for its own sake, at- 
tention has turned in recent years to the capabilities of 
LISA as a true astronomical observatory. Many papers 
[3l-[l9| have investigated just how well LISA can measure 
the parameters of the binaries it detects. This is often 
done using the Fisher-matrix method [13, HH , which es- 
sentially measures the local curvature of the posterior 
probability distribution for parameters in the vicinity of 
the maximum. Parameters for which the posterior is 
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more strongly curved (i.e., which more strongly affect 
the waveform) are measured more accurately than those 
for which the posterior is only weakly curved. Correla- 
tions between parameters are also extremely important. 
When two parameters are strongly correlated, it is difh- 
cult to "detangle" the influence of one on the waveform 
over the other. This means that the accuracy with which 
both parameters are measured is controlled by the one 
which is most poorly determined. 

Recent studies have considered how well parameters 
can be measured while doing the actual data analysis 
problem of removing signals from noise [22|. In effect, 
these studies simulate the parameter extraction process 
with enough detail to uncover issues such as multiple 
extrema of the posterior surfaces which are missed by 
the simpler (and cruder) Fisher analyses. Such "realis- 
tic" studies are typically much more CPU-intensive and 
cannot easily study parameter measurement issues over 
a broad swath of astrophysically important parameter 
space. Both families of studies have substantially ad- 
vanced our understanding of LISA's science reach in the 
past 5 or so years. 

Some parameters that are especially interesting are the 
intrinsic system properties, namely, the masses and spins 
of the black holes. Masses can be measured extremely 
well in the best cases [with a relative error of ^ 10~^ 
for individual masses and ~ 10^'"' for the chirp mass 
A4 = rr?^^ rr?2^ I (m\ -|- m'if-l^\. Spins are not measured 
quite so well but are still expected to be determined with 
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percent-level accuracy. By measuring these parameters 
for many systems, one can construct a merger history of 
black holes, and by extension, their host galaxies, learn- 
ing much about galaxy formation, black hole formation, 
AGN feedback, and so forth. 

We are also interested in measuring parameters extrin- 
sic to the system, namely, its position on the sky and 
its luminosity distance. With the position and distance 
(converted into an approximate redshift), astronomers 
can search the sky for probable electromagnetic counter- 
parts to the gravitational wave events. Various types of 
counterparts have been proposed, from signals during the 
inspiral '23'], to bright flashes at the time of merger [2^ 
(or even reductions in luminosity (25j). to long-delayed 
afterglows [2|| . The different scenarios arise because the 
behavior of gas around an inspiraling binary system is not 
well understood. A very different kind of electromagnetic 
counterpart can be produced by a kicked remnant black 
hole that triggers a telltale sequence of stellar disrup- 
tions [131. (Tidal disruption of stars may also allow us 
to flag the presence of a binary long before it enters the 
LISA band, allowing a better understanding of the space 
density of massive black hole binaries [Isj.) If a counter- 
part can be identified, the electromagnetic information 
can be combined with the gravitational information to 
reveal more about the astrophysics of the system. Coun- 
terparts may also make it possible for binary black holes 
to be used as probes of the cosmological distance-redshift 
relation, since the electromagnetic redshift and gravita- 
tional distance are determined independently [7|]. Unfor- 
tunately, finding a counterpart, even if a unique signature 
does exist, will not be easy. The typical error windows 
for LISA are ~ tens of arcminutes on a side at the end 
of inspiral, reduced from several square degrees in the 
weeks and months before merger. Still, this localization 
does give large survey telescopes like LSST (field of view 
~ 10 square degrees) a chance to study a particular area 
of the sky with advance warning [29 l| . 

One particularly important result in the study of 
LISA'S science capabilities was the discovery that includ- 
ing spin precession effects in the waveform model typi- 
cally improves the accuracy of parameter measurement 
I, a. Spin precession arises because of geodetic and 
gravitomagnetic general relativistic effects |30l |3T| . The 
orbital plane of the system also processes in order to pre- 
serve the total angular momentum on time scales shorter 
than the radiation reaction time. Together, these preces- 
sions modulate the amplitude and phase of the waveform, 
breaking correlations between certain sets of parameters 
and improving how well the members of those sets can 
be measured. The greatest improvement is to the mea- 
sured masses of a binary's members (accuracy typically 
improved by 1-2 orders of magnitude). The measured 
sky position angles and distance to the source are all im- 
proved by about half an order of magnitude, reducing the 
size of the sky position pixel in which one must search 
for a counterpart by a factor of ^ 10 (or the 3D voxel 
volume by a factor of ^ 30). 




The precession of one of the spins in a binary, to 1.5 
post-Newtonian order and averaged over an orbit, ^ is 
given by 



(1.1) 



where Si and S2 are the two spins, mi and m2 are the two 
masses, M = mi -I-TO2 is the total mass, /i = 11111712/ M is 
the reduced mass, L is the direction of the orbital angular 
momentum, and r is the orbital separation in harmonic 
coordinates. It is clear that precession is maximal when 
the spins of the system are orthogonal to the orbital angu- 
lar momentum and to each other and vanishes when the 
spins and orbital axis are aligned. In Q (hereafter Pa- 
per I), it was assumed that the relative orientation of the 
spins and the orbital angular momentum was completely 
arbitrary. The results of that paper are summarized in a 
series of histograms describing parameter measurement 
accuracy when the various angular momentum vectors 
are allowed to point in any direction. 

Recent studies have shown, however, that accreting gas 
in a system may evolve the spin in such a way that the 
spins are at least partially aligned with each other and 
with the orbit [H, [3^ . The degree of alignment depends 
on the temperature of the gas: In "hot gas" models, 
which have polytropic index 7 = 5/3, the spins align 
within 30° of the orbital axis. "Cold gas" models with 
7 — 7/5 align even more thoroughly, to within 10° [33]- 

Does spin-induced precession, now constrained by ini- 
tial conditions, still break degeneracies as efficiently as 
described in Paper I? Any degradation in parameter mea- 
surement capability could have a strong effect on the abil- 
ity to find electromagnetic counterparts. The results of 
M 01 13 be biased toward gas-free "dry" mergers, 
severely underestimating localization errors in gaseous 
"wet" mergers — the very systems which we are most 
likely to see electromagnetically. The effect of align- 
ment on mass and spin measurements is also interesting 
(though arguably less so, since even a factor of several 
degradation for these parameters would still imply excel- 
lent accuracy). 

The goal of this paper is to answer the question posed 
above. We do so with a Fisher-matrix analysis of pa- 
rameter measurement for binaries whose spins are par- 



^ The validity of averaging over an orbit can be quickly checked 
by comparing the precession time scale Tp^oc with the orbital 
time scale T^^h- For a-n equal-mass system, these time scales are 
roughly equivalent when r ~ 7M/8, where M is total mass. The 
orbit-averaged approximation is therefore quite good for most 
of the inspiral waveform but will begin to break down as we 
approach r ~ (2 — 3)M. At these radii, the post-Newtonian 
approximation breaks down as well, so our use of the precession 
equations there is best considered to be qualitatively indicative 
of the relevant physics, if not numerically exact. 
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tially aligned according to two wet merger models: hot 
gas, which aligns the spins and orbit to within 30°, and 
cold gas, which aligns to within 10°. We demonstrate 
that this degree of alignment can substantially degrade 
parameter accuracy but that one can "repair" much of 
this degradation by using the "full" waveform model, in- 
cluding harmonics beyond the leading quadrupole. In 
what follows, we will use the terms "gas-free" or "dry" 
interchangeably with the term "random spins," "hot gas" 
interchangeably with the phrase "spins aligned within 
30°," and "cold gas" interchangeably with the phrase 
"spins aligned within 10°." We also sometimes write "30° 
alignment" as shorthand for "spins aligned within 30°," 
and likewise for "10° alignment." (The two distributions 
contain systems with alignments less than 30° or 10°, al- 
though with a bias toward the upper end of the allowed 
range.) 

The outline of the paper is as follows. We begin in 
Sec. In] by describing the operation of our code, includ- 
ing the production of binary black hole waveforms, the 
LISA response, the noise model, and how we construct 
the Fisher matrix. We focus on changes from Paper I, 
leaving detailed description of the theory to that paper. 

In Sec. mil we then present results for parameter errors 
in wet mergers, examining both "hot" and "cold" models. 
We compare these results to the case of dry mergers, in 
which spin orientations are chosen to be completely ran- 
dom with respect to each other and to the orbital angular 
momentum. It should be emphasized that throughout 
Sec. Illli we consider only the leading quadrupole piece 
of the gravitational waveform (the so-called "restricted" 
post-Newtonian approximation). As expected, we find 
that spin alignment largely degrades LISA's ability to 
measure parameters. As a rough rule of thumb, we find 
that extrinsic parameters (sky position angles and lumi- 
nosity distance) are measured a factor of ^ 1.5 — 2 less 
accurately for 30° alignment and a factor of ^ 2 — 3 less 
accurately for 10° alignment. In the second case, align- 
ment eliminates most of the advantage gained by adding 
precession in Paper I. We find that the impact upon mea- 
sured masses and spins depends strongly on mass ratio, 
with degradation by a factor ^ 1 — 3 at 30° alignment, 
and a factor ^ 1 — 9 at 10° alignment. We find a hand- 
ful of cases in which partially aligned binaries actually 
do better than the randomly oriented systems. As we 
describe in Sec. IIIIl this is due to alignment increasing 
these systems' average SNR. 

To combat this degradation, we introduce another 
degeneracy-breaking effect. Much early work in LISA 
parameter estimation made use of the restricted wave- 
form model, in which only the quadrupole harmonic 
of the orbital phase was included and only the lead- 
ing "Newtonian" amplitude term was used with this 
harmonic (although the phase was constructed to high 
post-Newtonian order). This was done because the 
quadrupole harmonic dominates signal power, while the 
phase is the primary source of information about the 
signal. However, it has since been shown by several 



groups that including higher harmonics (and their post- 
Newtonian amplitudes, making the so-called full wave- 
form model) also breaks degeneracies and reduces pa- 
rameter errors [Tol-[T^. [T^ . The magnitude of the effect 
is comparable to the improvement seen by including spin 
precession. Recently, Klein et al. have presented an anal- 
ysis combining both spin precession and higher harmon- 
ics [16(1 . A similar analysis, based on an earlier version of 
our own code, was conducted by the LISA Science Team 
to investigate the science reach of the LISA mission; the 
results of this study are summarized in Ref. [s^] • 

In Sec. llVi we replace the leading quadrupole wave- 
form with the full waveform. For the case of random 
spins, our answers can be compared (with some caveats) 
to the results of [l^. We also compute the errors for 
wet, partially aligned binaries with the full waveform. 
When higher harmonics are included, parameter errors 
for partially aligned binaries are often no worse, or even 
better, than for the case of random spins and no higher 
harmonics. In these particular cases, higher harmonics 
can more than make up for the degraded impact of spin 
precession. We find this degree of improvement for the 
minor axis of the sky position error ellipse and for the 
luminosity distance in a majority of (mass) cases. The 
improvement is not quite so good for the major axis: Al- 
though higher harmonics can reduce errors by factors of 
^ 2 or more, this often does not completely make up for 
the loss of precession, especially at 10° alignment. Er- 
rors in the measured spin behave similarly to the major 
axis — their measurement is improved, but not enough 
to fully compensate for the impact of aligned spins. By 
contrast, we find that higher harmonics always improve 
mass measurements beyond what can be done with ran- 
dom spins alone. In fact, partial alignment in many cases 
improves mass measurements, thanks to increased SNR 
in these cases. 

We also briefly take a more detailed look at the rela- 
tive improvement from spin precession, higher harmonics, 
and their combination. We confirm previous results that, 
for extrinsic parameters, the impact of the combined ef- 
fects is not substantially greater than the impact of each 
effect alone. For mass errors, the higher harmonics dom- 
inate, with precession being almost irrelevant for the full 
waveform. For spin errors, however, the two effects do 
seem to be independent, with the combined improvement 
approximately equal to (or greater than) a simple multi- 
plication of the individual improvements. 

We conclude in Sec.|V]by summarizing our results and 
discussing additional studies that must be done before 
the question of LISA parameter estimation is fully un- 
derstood. Throughout this paper, we use geometrized 
units in which G = c — 1. A useful conversion factor is 
that 1 M© = 4.91 X 10"*^ s = 1.47 km. 

Since this paper was originally written, budget con- 
straints have caused a rescoping of the LISA mission, 
and the mission that eventually flies may differ from the 
"classic" configuration considered here. We continue to 
focus our analysis on measurements using LISA Classic 
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for two reasons. First, the design of the rescoped mis- 
sion is in flux. Until a design is fixed and its associated 
sensitivity known, we cannot study how weh it wih make 
measurements. Second, our goal is to make comparisons 
with previous studies that were based on the classic de- 
sign. As such, it is most appropriate for us to use this 
design as well. We note that our conclusions should be 
robust in the sense that the general trends we find re- 
garding the impact of spins and higher harmonics will be 
relevant to any LISA-like design (at least for designs that 
have five or six links, so both waveform polarizations can 
be simultaneously measured). 



II. PARAMETER ESTIMATION CODE 

The code used in this paper is a version of the 
MONTANA-MIT codc uscd by the LISA Parameter Esti- 
mation Taskforce [3^ , updated with some new features 
and bug fixes. In this section, we describe the relevant 
features of the code, especially how it differs from the 
code of Paper I We refer the reader to Paper I for 
more detailed discussion of the waveform and parameter 
estimation theory. 



A. Massive black hole binary waveform 

The waveform from a massive black hole binary coales- 
cence is traditionally divided into three distinct phases: 

(1) the inspiral of the two holes, which can be described 
by the post-Newtonian expansion of general relativity; 

(2) the merger of the two holes into a common event 
horizon, describable only by full numerical relativistic 
simulations; and (3) the ringdown of the final hole into 
the stationary Kerr solution, which can be described by 
black hole perturbation theory. In this work, we con- 
sider only the inspiral, which for LISA sources can last 
for months to years, accumulating large amounts of SNR 
and parameter information. Because of this fact, as well 
as the ease of using the post-Newtonian approximation, 
inspiral-only waveforms have traditionally been used in 
most, though not all, LISA parameter estimation stud- 
ies. Ringdown information was first studied on its own by 
Berti, Cardoso, and Will More recently, Mc Williams 
et al. added both the merger and ringdown to the in- 
spiral, albeit for nonspinning binaries with an a priori 
known mass ratio [13, Il3 ■ They showed that the merger 
can add a significant amount of parameter information, 
about a factor of 3 improvement in measurement accu- 
racy for all parameters but mass. Work in progress will 



^ We note that the codes used in |34| were all found to produce 
the same answers provided they used the same noise models, the 
same signal cutoffs, and so on. Our code has been well-tested 
in as much as other codes had the same features to compare 
against. 



consider the impact of an unknown mass ratio, as well as 
spins. 

The inspiral waveform can be described by 17 param- 
eters: the masses of the black holes, mi and 7712; their 
dimensionless spins, xi = ISil/mf and X2 = |S2|/'Ti2; 
the spin angles at some particular reference time to, 
cos 6*51 (io), cos 6*52 (to), (t^sAto), and ^52(^0); the orien- 
tation angles of the orbital angular momentum at to, 
cos0L(io) and 0L(to); the eccentricity e; the periastron 
angle 7; the position of the binary on the sky, cosO^ 
and 0jv; the luminosity distance D^; a reference time 
trof (possibly different from to); and a reference phase 
$ref = *i'(trof)- In this work, we assume quasicircular 
orbits, eliminating e and 7 and reducing the parameter 
set to 15. This assumption is also quite common, since 
radiation reaction has long been expected to circularize 
binaries [sH]. It should be noted, however, that recent 
studies indicate that gas [H, [s^l and/or steUar interac- 
tions [3^ may cause binaries to retain a small, but sig- 
nificant, residual eccentricity when they enter the LISA 
band. Recent work by Key and Cornish investigates 
the impact of this residual eccentricity using a nontrivial 
extension of our code. 

In Paper I, we used the post-Newtonian parameters tc 
and (j)c as the reference time tj-of and phase 'I'lef- These 
parameters are, respectively, the time and phase when 
the post-Newtonian frequency formally diverges. How- 
ever, Paper I made a slight error in determining the post- 
Newtonian frequency and phase. To understand this er- 
ror and how to correct it, begin with the time derivative 
of orbital angular frequency w = 27r/orb (shown here to 
second post-Newtonian order) 
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where 77 = pi/M is the reduced mass ratio, /3 is a spin- 
orbit coupling term, and cr is a spin-spin coupling term. 
Exact expressions for /3 and a are given in Paper I. Equa- 
tion (12. ip must be integrated once to obtain Lo{t) and 
twice for the orbital phase $orb(t)- When the spins do 
not precess, this integration can be done analytically to 
some specified post-Newtonian order. In Paper I, the 
analytic results were used, but with the time-dependent 
expressions for /3 and a plugged in at the end of the pro- 
cess. This is technically not correct: The time-dependent 
spins should be inserted into (|2.ip , and then that expres- 
sion should be numerically integrated to produce uj(t) 
and '^orh(t). This is not difficult, only requiring two ad- 
ditional differential equations in the Runge-Kutta solver 
of Paper I. However, it means that tc and $c are no longer 
acceptable references, since the numerical integrator can- 
not reach infinite frequency. We describe our current ap- 
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proach momentarily. 

Another change from the code used in Paper I is in 
the choice of cutoff frequency for the inspiraL In Pa- 
per I, the inspiral was stopped at the frequency of the 
Schwarzschild innermost stable circular orbit (ISCO), 
r = 6M. This assumption is poor for two reasons. First, 
while r = 6M is the ISCO for a test particle orbiting a 
single Schwarzschild hole of mass A/, the dynamics of the 
two-hole system are much more complex, and the transi- 
tion to plunge and merger is not so well-defined. Second, 
we are considering Kerr black holes, for which even in the 
point-particle limit the innermost stable orbit can vary 
from r = 9M to r = M depending on the spin of the 
hole, with a concomitantly wide variation in the ISCO 
frequency. A better solution is to stop the inspiral at the 
minimum energy circular orbit (MECO), the orbit which 
minimizes the expression for post-Newtonian energy [s^ : 



3M2 



1 + |^)L.S, 

4 mi 
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4 7712 
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(Si •S2-3(L-Si)(L-S2)) 



(A/c.)4/3 
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The MECO is known to be a better approximation to the 
inspiral-plunge transition than the ISCO, and it properly 
takes spins into account. 

Using the MECO gives us a better reference point for 
our time and phase than the coalescence time and phase 
tc and described above. We choose to — trcf — iMECO 
and $rcf = 'I'meco and then integrate the spin, fre- 
quency, and phase evolution equations backwards from 
the MECO lot = Q. The backwards integration provides 
stability in the Fisher-matrix calculation: We align the 
waveforms when they are largest, thus making it easier 
to introduce slight perturbations. 

As seen in Eq. (12. ip . we calculate the phase out to sec- 
ond post-Newtonian (2PN) order. (By numerically inte- 
grating ()2.1|) to obtain ijo{t) and $orb(0: specifically 
are choosing the "TaylorT4" PN approximant [1^.) We 
integrate the spin precession equations out to 1.5PN or- 
der, which includes IPN spin-orbit and 1.5PN spin-spin 
terms. It is worth noting that all of the relevant quanti- 
ties are known to higher post-Newtonian order. Work in 
preparation shows that including terms beyond the order 
we include here only causes a slight quantitative change 
in the accuracy with which parameters are measured [4l| . 
In Sec, mil we use the restricted post-Newtonian approx- 
imation, in which we only consider the quadrupole term 
($ = 2$orb) with its lowest order, Newtonian amplitude. 
In Sec. IIV| we use the full post-Newtonian waveform, 
which includes all harmonics to 2PN order in amplitude. 



B. LISA response and noise 



The LISA response used in this paper differs from the 
response used in Paper I. For signals which do not reach 
above 3 X 10^^ Hz, we use the same low-frequency approx- 
imation used in that paper. In this approximation, we 
ignore the transfer functions which arise due to the finite 
arm lengths of the detector. This approximation is very 
inaccurate above / ~ 3 x 10~^ Hz. With the addition of 
higher harmonics, many signals now reach into this range 
where the transfer functions become important. 

The full LISA detector response is somewhat compli- 
cated to model. Three existing codes provide the full 
response: the LISA Simulator [l^]. Synthetic LISA ^1], 
and LIS ACode [43| . Interfacing with one of these codes 
would significantly slow our analysis, making it difficult 
to perform large Monte Carlo studies over our param- 
eter space. We seek a simpler response function which 
includes the finite arm length transfer functions but ig- 
nores some of the more complicated issues. 

The three LISA spacecraft follow eccentric orbits 
around the Sun at 1 AU. The individual orbits combine 
in such a way that the LISA constellation maintains, at 
first order in orbital eccentricity, an equilateral triangle 
formation. By going beyond this leading order, one finds 
that the arm lengths vary by a small amount on month- 
long time scales. The variation in LISA arm lengths is 
the reason for the development of time delay interfer- 
ometry (TDI) techniques [i^ to eliminate laser phase 
noise, which cancels exactly in equal-arm interferometers 
like LIGO. In our detector model, we approximate the 
constellation as having arm lengths that are equal at all 
times. Our model detector is a "rigid" equilateral trian- 
gle. 

The other complexity in the full LISA response is that 
the spacecraft move during the measurement, causing 
"point-ahead" effects which must be taken into account. 
We assume instead an "adiabatic" detector, in which for 
each time that we require the detector response, the de- 
tector is considered to be motionless for that time. The 
spacecraft then adiabatically move to their next position 
for the next sample point. This rigid, adiabatic approxi- 
mation is known to be equivalent to the full response up 
to very high frequency (~ 500 mHz) and thus is appro- 
priate for our Fisher-matrix analysis p^ . 

For the rigid, adiabatic approximation, the code pro- 
duces Michelson variables X,Y , and Z as defined in [43], 
Eqs. (10) and (11). Note that these are technically not 
TDI variables. Since we do not have to subtract phase 
noise, there is no need to include another pass through 
the interferometer (cf. the "real" TDI variables in Eq. 
(13) of [4^). They do contain the same information, 
though, so we may refer to them as (pseudo, equal-arm) 
TDI Michelson variables in this paper. From them, we 
can construct noise-orthogonal TDI variables A, E, and 
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T, defined as 



-{2X - Y - Z) , 
3 

T =\{X + Y + Z). 



(2.3) 
(2.4) 

(2.5) 



A, E, and T are used to calculate the SNR and the Fisher 
matrix. Note that, as defined in Eqs. (10) and (11) of [43j. 
these are fractional-frequency variables. We can convert 
them to equivalent strain by integrating the signal in the 
frequency domain and then multiplying by c/(47rL). For 
the low-frequency case, we use the Michelson signal X 
and the noise-orthogonal signal {X + 2Y)/y/3. These 
combinations are denoted hj and hn in Paper I (which 
in turn follows the convention of Cutler Q). The low- 
frequency approximation is constructed so that these sig- 
nals are already expressed as equivalent strain. 

The LISA noise power spectral density S'„(/) com- 
prises two parts, instrumental noise and confusion noise 
due to unresolved white dwarf binaries in the Galaxy. 
Instrumental noise consists of both position noise, due 
to photon shot noise and other effects along the optical 
path, and acceleration noise, due to proof mass motion. 
The total instrument noise in the A and E (strain) chan- 
nels is given by 
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and the T (strain) noise is given by 
1 



3L2 
1 



(1 - cosx) Sp{f) 
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Here L = 5x10^ km is the LISA arm length, x = 27r/L / c, 
Sp(f) = 3.24 X 10~22 m^/Hz is the position noise bud- 
get, and Saif) — 9 X 10^^° m^/s^/Hz is the accelera- 
tion noise budget. For the low- frequency approximation, 
we can calculate similar expressions for the two noise- 
orthogonal channels and then take cos a; — 1, although 
this attention to detail makes little difference for the fre- 
quencies of interest. Notice that the position noise and 
acceleration noise are both assumed to be white, with no 
frequency dependence. However, because it is expected 
that LISA'S acceleration noise performance will degrade 
somewhat from this white form below 10^^ Hz, we have 
also added a "pink" acceleration noise term, with a slope 
oif-'. 

Confusion noise is constructed from the residuals of a 
fit to the Galaxy i46i] in the Mock LISA Data Challenge 



[23 [. An approximate analytic expression for the confu- 
sion noise can be found in JJ,], Eq. (10). It is added to 
instrument noise for the A and E channels (or the orthog- 
onal low-frequency channels) to obtain the total noise. It 
is not added to the T channel because it occurs only at 
low frequency, where that channel adds nothing to the 
analysis. 

Finally, although it is not expressed explicitly in (|2.6p 
and (|2.7p . we enforce a low-frequency cutoff of 3 x 10~''^ 
Hz and do not include any contribution from the signal 
below that frequency in our analysis. (This frequency is 
the lowest frequency at which LISA is planned to have 
good sensitivity to gravitational waves; though it will 
have sensitivity to sources at lower frequencies, the noise 
characteristics below / = 3 x 10~^ Hz cannot be guaran- 
teed.) 



Construction of the Fisher matrix 



The Fisher matrix Fi, is defined as 



dh 



dh 



where h is the gravitational wave signal, 
parameters which describe it, and 



(2.8) 



are the 15 



(a|6) =4Rc / df 



r a*if)bif) 

Snif) 



(2.9) 



is a noise-weighted inner product. The inverse of the 
Fisher matrix is the covariance matrix, which contains 
squared parameter errors along the diagonal and correla- 
tions elsewhere. To calculate the Fisher matrix, we need 
the waveforms in the frequency domain. In Paper I, we 
actually did all calculations in the frequency domain by 
using the stationary phase approximation. This approxi- 
mation relies on a separation of time scales and is known 
to be quite good for nonspinning binaries, where the in- 
spiral time scale Tjnsp is much larger than the orbital 
time scale Toib- However, when precession is included in 
the waveform, an additional time scale Tproc comes into 
play, with Ti„sp > Tprcc > Toib- We have seen that with 
precession, the stationary phase approximation tends to 
smooth out sharp features in the Fourier transform, po- 
tentially reducing the information content. The problem 
becomes worse as the impact of precession increases (i.e., 
with higher spin values, or for highly nonaligned spins 
and orbit). To avoid introducing any errors due to this 
approximation, we here calculate our waveforms in the 
time domain and then perform a fast Fourier transform 
(FFT) to bring them into the frequency domain. 

This approach has two major limitations. First, it is 
much slower than the stationary phase approach, since we 
need to calculate many time samples to observe Nyquist 
sampling requirements and we then need to compute the 
FFT. Second, the FFT assumes a periodic signal. Be- 
cause we have a finite signal which looks much different 



7 



at the end than at the beginning, we must introduce some 
kind of window in order to taper the signal to zero at the 
beginning and end. We use a Hann window (actually half 
a Hann window at each end of the signal) . This window 
substantially reduces "ringing," or spectral leakage prob- 
lems. However, it also cuts out part of the signal. This is 
particularly unfortunate for the strongly chirping inspi- 
ral, since much of the signal power is contained in the last 
few cycles. By windowing the signal, we lose some of this 
power. This may cause our SNR and errors to be smaller 
and larger, respectively, than they would be for a "real," 
physical signal. The best solution to this problem would 
be to include the merger and ringdown portions of the 
signal, allowing it to fade to zero in a physical, not artifi- 
cial, way. For now, we must simply accept the windowing 
as part of the definition of the (unphysical) inspiral-only 
waveform. 



III. PARAMETER ESTIMATION IN 
PARTIALLY ALIGNED BINARIES: ONLY THE 
QUADRUPOLE HARMONIC 

Here we describe the parameter estimation capabili- 
ties of LISA without including the influence of higher 
harmonics. In order to consider a wide range of LISA 
sources, we choose only three parameters explicitly, the 
two masses of the system and the luminosity distance. 
For the masses, we consider a variety of systems ranging 
in total mass from 2 x 10^ Mq to 2 x lO'^ Mq, with a mass 
ratio from 1 — 10. On the other hand, we consider only 
sources at z = 1, corresponding to a luminosity distance 
of 6.64 Gpc using our choice of cosmological parameters. 
Errors at other redshifts can be constructed using the 
errors at z = 1. We note that the results at masses mi 
and 1712 and redshift z' can be simply related to the re- 
sults at masses mi{l + z')/{l + z) and m2(l + z')/{l + z) 
and redshift z. This is because all time scales in the sys- 
tem are derived from the masses. Since time scales are 
lengthened (frequencies shortened) by the cosmological 
redshift, a binary at higher redshift behaves like a binary 
at lower redshift but with a higher mass. The quantity 
to(1 -f z) is generally called the redshifted mass, where 
m is the mass measured locally at the rest frame of the 
binary. (When we quote masses in this paper, we always 
mean the rest-frame mass, remembering that when put 
into the waveform formulas of Paper I, they must be mul- 
tiplied by 1 + z.) The amplitude of the waves at redshift 
z' is decreased by a factor ^ = Dl{z)/Dl{z') over the 
corresponding binary (i.e., the binary with the same red- 
shifted mass) at redshift z. This increases the errors by 
1/^ over that corresponding signal. 

The other 12 parameters of the system are generated 
essentially at random, with 1000 different Monte Carlo 
realizations. For example, iiviECO is chosen from within 
an assumed three- year mission time, meaning that some 
early binaries will have abnormally short signals for a 
given mass. Spin magnitudes are chosen uniformly be- 



tween and 1, and <I>meco is chosen uniformly between 
and 27r. Cosines of angles are chosen uniformly between 
— 1 and 1, while longitudinal angles are chosen uniformly 
between and 27r. In the case of random spins (as in 
Paper I), the procedure is then complete. In the case of 
partially aligned spins, the main focus of this paper, we 
use the randomly generated parameters to integrate the 
spin precession equations backwards from the MECO to 
t = 0. We assume that any alignment at t = is solely 
due to gas. If either of the resulting spin-orbit angles is 
greater than the model's restriction (30° or 10°), we ran- 
domly select new spin orientation angles (at MECO) and 
try again. This procedure guarantees that all sources will 
have the desired amount of alignment at the start of the 
signal. However, our sample will include some sources 
(~ 30%) which move out of alignment by MECO. Since 
these sources generally precess more strongly than the 
others, they tend to improve the overall distribution of 
parameter errors, especially for spin magnitude. 

Figure [T] shows a histogram of the Monte Carlo results 
for a binary with mi = lO^M© and TO2 = 3 x 10^ Mq. 
We show the major axis of the sky position error ellipse, 
2a, comparing the cases of randomly aligned spins to 
spins restricted to be aligned within 30° (for hot gas) 
and 10° (for cold gas) of the orbital angular momentum. 
We see that partial alignment of the spins and orbital 
angular momentum degrades EISA's localization capa- 
bility. For the partially aligned cases, the shapes of the 
histograms resemble the strongly peaked "no precession" 
results of Paper I more than the roughly flat random- 
spin histogram. The medians of the distributions also in- 
crease: While randomly oriented binaries have a median 
major axis of 34.8 arcminutes, systems aligned within 30° 
have a median 2a of 62.3 arcminutes. For 10° alignment, 



140 




2a (arcminutes) 



FIG. 1: Distribution of 2a, the major axis of tlie sky posi- 
tion error ellipse, for binaries with randomly aligned spins 
(dotted line), spins restricted to within 30° of the orbital 
angular momentum (dashed line), and spins restricted to 
within 10° of the orbital angular momentum (solid line). Here 
mi = 10® Mq, m2 = 3 X IO^Mq, and z = 1. 
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FIG. 2: Same as Fig. [Tl but for the minor axis 2b. 



this degrades further to 90.5 arcminutes. This is a factor 
of 2.6 degradation from the case of random aUgnment, 
just short of the factor ~ 3 improvement seen in Paper I 
when precession is introduced into the waveform model. 
In essence, by restricting the spin angles to within 10° 
of the orbital angular momentum, we have eliminated 
almost all of the advantage gained from including pre- 
cession effects in the waveform. 

Figure [2] shows results for the minor axis of the sky 
position error ellipse,'* 2b. The results are similar: The 
median value of 2& increases from 24.6 arcminutes for 
random spins to 40.6 arcminutes for 30° alignment and 
to 58.6 arcminutes for 10° alignment. Together with the 
results for the major axis, these numbers imply that the 
total sky position area increases by a factor > 6 when 
binaries have closely aligned angular momentum vectors, 
strongly impacting the ability of LISA to find electro- 
magnetic counterparts to the GW signal. 

Table U shows the major and minor sky position axes 
for a range of masses, in the random-spin, hot gas, and 
cold gas cases. We see that degradation of a factor ~ 2—3 
between the "no gas" and "cold gas" (10° alignment) 
cases occurs rather consistently for different masses and 
mass ratios. 

The other extrinsic parameter of interest is the lumi- 
nosity distance D^. Figure [3] shows the fractional errors 
in Dl for different degrees of spin alignment. Again, we 
see that restricting the spin angles dramatically affects 
measurement: The median of 5.24 x lO^'^ for random 
spin orientation doubles to 1.01x10^^ when the spins are 
aligned within 30° and nearly triples to 1.36 x 10~^ when 
the spins are aligned within 10°. However, this particu- 
lar degradation is almost certainly immaterial, at least at 



mi (Mq) 


m2(M0) 


No 
2a 


gas 
2b 


Hot 
2a 


gas 
2b 


Colo 
2a 


gas 
2b 


10'^ 


10^ 


27.0 


16.4 


40.7 


25.7 


53.8 


35.1 


3 x 10^ 


10^ 


17.5 


11.7 


30.1 


17.9 


53.7 


34.8 


3 X 10^ 


3 X 10^ 


33.3 


19.0 


45.9 


27.1 


63.1 


42.4 


10'' 


10^ 


23.3 


18.3 


35.9 


21.6 


61.6 


38.4 


lO'* 


3 X 10^ 


34.8 


24.6 


62.3 


40.6 


90.5 


58.6 


10^ 


10® 


56.9 


37.5 


87.7 


57.2 


105 


68.3 


3 X 10® 


3 X 10^ 


39.0 


33.6 


57.0 


36.8 


105 


68.1 


3 X 10® 


10® 


45.5 


32.4 


83.3 


49.0 


131 


77.9 


3 X 10® 


3 X 10® 


71.9 


43.6 


126 


75.6 


168 


106 


10^ 


10® 


47.3 


40.2 


70.7 


46.6 


132 


83.7 


10^ 


3 X 10® 


67.3 


45.3 


131 


75.7 


234 


143 


10^ 


10^ 


160 


84.8 


281 


136 


581 


323 



TABLE I: Median sky position major axis 2a and minor axis 
2b, in arcminutes, for binaries of various masses at z — 1, in 
the "no gas" (random-spin), "hot gas" (30° alignment), and 
"cold gas" (10° alignment) cases. 



low redshift, since the error remains much smaller than 
the ^ 5% error produced by weak gravitational lensing 
at z ~ 1. For sources at higher redshift, this degrada- 
tion may be more important. Table HI] shows luminosity 
distance errors for different masses. Like the sky posi- 
tion, the degradation is about a factor of ^ 1.5 — 2 for 
30° alignment and ~ 2 — 3 for 10° alignment. Note that 
for the larger masses we consider, the degradation pushes 
the GW distance error to a value comparable to or even 
larger than the weak lensing error. 

We now turn to the intrinsic parameters of the sys- 
tem, its masses and spins. Figs. |4] and [5] show the 
errors in the two black hole masses for the three cases 
we consider. Medians of Ami /mi are 4.84 x 10^'' for 
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10 

ADl/Dl 



* Note that unlike in Paper I, the results for 2b do not feature a 
long tail of small errors. We have confirmed that this effect was 
caused by a bug in the code used in Paper I. 



FIG. 3; Same as Fig. [T] but for the fractional error in lumi- 
nosity distance, ADl/Dl. 
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mi (Mq) 


m2(M0) 


No gas 


Hot gas 


Cold gas 


10^ 


10^ 


4.16 X 10^2 


7.97 X IQ-^ 


0.0130 


3 X 10^ 


10^ 


2.59 X 10"^ 


6.03 X IQ-^ 


0.0101 


3 X lO'^ 


3 X 10^ 


5.54 X 10~^ 


9.23 X 10"^ 


0.0121 




10^ 


3.67 X 10"^ 


5.92 X IQ-^ 


0.0113 


10" 


3 X lO'^ 


5.24 X 10"^ 


0.0101 


0.0136 


10" 


10" 


9.37 X 10"^ 


0.0137 


0.0175 


3 X 10" 


3 X 10^ 


5.58 X 10"^ 


9.19 X 10"^ 


0.0147 


3 X 10" 


10" 


7.06 X 10"^ 


0.0130 


0.0191 


3 X 10" 


3 X 10" 


0.0137 


0.0207 


0.0279 


10^ 


10" 


7.67 X 10"^ 


0.0135 


0.0242 


10^ 


3 X 10" 


0.0129 


0.0243 


0.0429 


10^ 


10^ 


0.0441 


0.0613 


0.0974 



TABLE II: Same as Table U but for the fractional error in 
luminosity distance, ADl/Dl- 
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FIG. 5: Same as Fig. [T] but for the fractional error in mass, 
Am2/m2. 



the random-spin case, 5.70 x 10^'^ for a system with hot 
gas, and 8.23 x 10~^ for a system with cold gas. For 
Am2/m2, these numbers are 3.84 x 10~^, 4.54 x 10^'^, 
and 6.55 x 10"'^, respectively. The impact of partially 
aligned spins does not seem to be as strong on the masses 
as on the sky position; the mass errors change by less 
than a factor of 2. In Paper I, we looked at preces- 
sion improvements not in individual masses but in chirp 
mass and reduced mass, where we saw factors of ^ 10 
and ~ 100 — 1000 improvement, respectively. Clearly, 
restricting the spin directions does not remove this en- 
tire improvement; even a limited amount of precession 
appears to significantly aid mass determination. We can 
check this assertion using our new code by running a case 
with only 1° alignment between the spins and the orbit. 
We find that the 10° results improve on the 1° results by 
a factor of 3. By contrast, the sky position and distance 
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FIG. 4: Same as Fig. [T] but for the fractional error in mass. 
Ami /mi. 



errors differ by only 15 — 20%. 

Table UlIl shows the results for different masses. We see 
a much stronger dependence on mass ratio here than for 
the extrinsic parameters. For example, while the cold gas 
degradation is less than a factor of 2 for the (roughly) 3:1 
mass ratio case considered in Figs. H] and [SJ it reaches a 
factor of ^ 9 for the equal- mass case mi = TO2 — 10^ Mq. 
This is unusual, since precession is known, at least for ex- 
trinsic parameters, to have a stronger impact for unequal 
masses due to increased complexity in the signal. It is 
possible that the lack of this complexity essentially "gives 
away" that the masses are equal, making them easier to 
determine from the extremely well-measured chirp mass. 

Interestingly, there are some examples of 10:1 mass ra- 
tio systems that break our general trend; in these cases, 
we find that partially aligned spins actually do better 
than random spins. This seemingly counterintuitive re- 
sult can be explained by our choice of the minimum en- 
ergy circular orbit (MECO) as the waveform cutoff. Bi- 
naries with aligned spins have a smaller MECO (with 
a corresponding high inspiral cutoff frequency) and thus 
accumulate more SNR than those with spins out of align- 
ment (as many in the random-spin sample will be). Fig- 
ure [5] shows the SNR for all three cases in a 10:1 binary. 
We see that the SNR is substantially larger for the par- 
tially aligned cases (medians of 2588 and 2592, for 30° 
and 10°, respectively) than the randomly aligned case 
(median of 1445). Even though these binaries precess 
less, the increase in SNR makes up for it in parameter 
estimation. It is worth noting that this effect could also 
be of use in detecting and measuring particularly high- 
mass binaries. For randomly chosen spins, such a binary 
might be mostly or completely out of the LISA band. 
However, if the spins are aligned by interactions with 
gas, the MECO frequency will be pushed into band. 

Finally, we consider how restriction of spin angles af- 
fects measurement of spin magnitudes; Figs.[7]and[S]show 



10 



nil (Mq) 


m2{M@) 


No 
Ami/mi 


gas 
Am2/m2 


Hot 

Ami /mi 


gas 
Am2/m2 


Cole 
Ami /mi 


gas 

Am2/m2 


10^ 


10^ 


3.23 X 10"^ 


3.24 X 10"^ 


9.84 X 10"=^ 


9.84 X 10"^ 


0.0284 


0.0284 


3 X 10^ 


10^ 


3.02 X 10"^ 


2.45 X 10"^ 


3.95 X 10-=^ 


3.20 X 10"^ 


5.94 X 10"^ 


4.83 X 10"^ 


3 X 10^ 


3 X 10^ 


4.50 X 10"^ 


4.50 X 10"^ 


0.0129 


0.0128 


0.0327 


0.0328 


10^ 


10^ 


2.72 X 10~^ 


1.90 X 10"^ 


1.90 X 10-=* 


1.32 X 10"^ 


2.22 X 10"^ 


1.55 X 10"^ 


lO'' 


3 X 10^ 


4.84 X 10"^ 


3.84 X 10^^ 


5.70 X 10"=* 


4.54 X 10~^ 


8.23 X 10"^ 


6.55 X 10'^ 


10'^ 


10® 


8.05 X 10"^ 


8.05 X 10^^ 


0.0197 


0.0197 


0.0475 


0.0475 


3 X 10® 


3 X 10^ 


5.81 X 10"^ 


4.01 X 10"^ 


4.49 X 10"^ 


3.07 X 10"^ 


5.10 X 10"^ 


3.51 X 10"^ 


3 X 10® 


10® 


0.0121 


9.73 X 10^^ 


0.0165 


0.0132 


0.0233 


0.0187 


3 X 10® 


3 X 10® 


0.0239 


0.0237 


0.0536 


0.0533 


0.109 


0.109 


10^ 


10® 


0.0176 


0.0118 


0.0167 


0.0109 


0.0205 


0.0133 


10^ 


3 X 10® 


0.0431 


0.0336 


0.0581 


0.0447 


0.0924 


0.0713 


10^ 


10^ 


0.381 


0.388 


0.424 


0.423 


0.967 


0.971 



TABLE III: Same as Table H] but for the mass errors Ami /mi and Am2/m2. 
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FIG. 6: Signal-to-noise ratio for binaries with randomly 
aligned spins (dotted line), spins restricted to within 30° of 
the orbital angular momentum (dashed line), and spins re- 
stricted to within 10° of the orbital angular momentum (solid 
line). Here mi = 10® M©, m2 = IO^Mq, and 2 = 1. 




FIG. 7: Same as Fig. \T\ but for the error in spin magnitude, 
Xi = |Si|/m?. 



FIG. 8: Same as Fig. [T] but for the error in spin magnitude, 
X2 = |S2|/mi. 



these results. For xi, the median varies from 4.55 x 10""^ 
for no gas to 5.38 x lO"'^ for hot gas to 1.31 x 10~^ for 
cold gas. For X2, the situation is similar; the medians are 
1.48x10-2 (no gas), 1.72x10-2 (hot gas), and 5.11 x lO'^ 
(cold gas). Although the spin errors are degraded by par- 
tial ahgnment, the amount of degradation is somewhat 
curbed by the contribution of binaries which precess away 
from alignment before MECO. In addition, the errors at 
10° alignment are roughly an order of magnitude bet- 
ter than at 1° alignment. Similar to the situation with 
mass measurements, even a small amount of precession 
can have a huge impact on measuring spin. 

Table IIVI gives spin errors for a broader range of 
masses. Like the mass errors, there is a strong depen- 
dence on mass ratio. In this case, however, the worst 
degradation occurs not for equal masses, but for 3:1 mass 
ratios, with factors of up to 4 increases in Axi and factors 
of up to 6 increases in Ax2- The 10:1 cases show some 
degradation at 10° alignment, but many cases are slightly 
improved at 30° alignment. As in the case of mass mea- 
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mi (A/q) 


m2(M0) 


No 

Axi 


gas 

AX2 


Hot ga 
Axi 


Ax2 


Cold g 
Axi 


as 

AX2 


10^ 


10^ 


0.0217 


0.0210 


0.0311 


0.0310 


0.0391 


0.0388 


3 X 10^ 


10^ 


3.18 X 10"^ 


7.61 X 10"^ 


4.98 X 10"^ 


0.0155 


0.0125 


0.0467 


3 X 10^ 


3 X 10^ 


0.0321 


0.0315 


0.0390 


0.0382 


0.0430 


0.0417 


10^ 


10^ 


1.30 X 10"^ 


0.0355 


1.05 X 10"^ 


0.0225 


2.22 X 10"^ 


0.0358 


10^ 


3 X 10^ 


4.55 X 10"^ 


0.0148 


5.38 X 10"^ 


0.0172 


0.0131 


0.0511 


lO'' 


10« 


0.0534 


0.0505 


0.0655 


0.0645 


0.0774 


0.0753 


3 X 10*^ 


3 X 10^ 


2.38 X 10"^ 


0.0601 


1.70 X 10"^ 


0.0499 


3.24 X 10"^ 


0.0659 


3 X lO'* 


10« 


9.92 X 10"^ 


0.0186 


0.0111 


0.0264 


0.0252 


0.0787 


3 X 10^ 


3 X 10" 


0.134 


0.125 


0.197 


0.198 


0.234 


0.232 


10" 


10" 


4.86 X 10"^ 


0.124 


3.69 X 10^^ 


0.180 


7.03 X 10-^ 


0.231 


10" 


3 X lO" 


0.0266 


0.0446 


0.0298 


0.0852 


0.0652 


0.191 


10" 


10" 


1.69 


1.54 


1.28 


1.31 


1.86 


1.87 



TABLE IV: Same as Table m but for the spin magnitude errors Axi and Ax2. 



surements (cf. Table Ullj) . this can be attributed to the 
increased SNR for aligned binaries. 



IV. PARAMETER ESTIMATION IN 
PARTIALLY ALIGNED BINARIES: INCLUDING 
HIGHER HARMONICS 

We now move beyond the leading quadrupole wave- 
form to the full waveforms. Post-Newtonian correc- 
tions to the waveform amplitude are included up to 
2PN order, including both additional quadrupole terms 
((j) = 2<i>orb) and subleading ("higher") harmonics be- 
yond the quadrupole. For example, the barycentric wave- 
form h^(t) can be written up to IPN order in the ampli- 
tude as 



h^{t) =^|(1 + cj) cos2<I>o,, + x^/^|^[(5 + c^) 



X cos $orb - 9(1 + C^) COS 3$orb] + X 



(19 



9cf - 2cf - r;(19 - llct - 6cf)) cos2$orb 
4 



-5^(1 + c2)(l - 377) C0s4$orb 



(4.1) 



where x 



[Mu 



^2/3 



cos(L • n), Si 



sin(L 



and Am — mi — 777,2. Here we see both extra harmonics 
and a IPN correction to the quadrupole harmonic. Note 
that the odd harmonics only contribute if mi ^ 771,2; just 
like spin precession, higher harmonic corrections are more 
complex for unequal masses. Further terms (including 
the X polarization) can be found in (albeit with some 
differences in sign convention). 

It has been shown that higher harmonic corrections can 
improve parameter estimation much like spin precession 



does [3 • In the case of higher harmonics, degen- 

eracies are broken due to the different sky position depen- 
dence of each harmonic. However, these studies did not 
include precession and so could not comment on how the 
two effects would combine. More recently, both effects 
were included in a parameter estimation study by Klein 
et al. (Ref. [l^). Their results demonstrate that includ- 
ing both precession and higher harmonics improves mea- 
surement accuracy, but, at least for extrinsic variables 
(sky position and distance), the combined improvement 
is not as drastic as the improvement from each effect on 
its own. This indicates that at least in some cases, pre- 



80 
70 
I 60 
a 50 



10 




10" 



10' 10' 
2a (arcminutes) 



10- 



10' 



FIG. 9: Distribution of 2a, the major axis of the sky posi- 
tion error ellipse, for binaries with randomly aligned spins 
(dotted line), spins restricted to within 30° of the orbital 
angular momentum (dashed line), and spins restricted to 
within 10° of the orbital angular momentum (solid line). Here 
nil ~ 10" Mq, 1712 = 3x IO^A/q, and z = 1. Higher harmonics 
are now included in the waveform model. 



cession and higher harmonics encode similar information. 
We might therefore expect that in partially ahgned bina- 
ries for which spin precession exists but is suppressed, 
the inclusion of higher harmonics may make up for this 
suppression and restore much, if not all, of the lost pa- 
rameter estimation capability. In this section, we test 
that expectation. 

Figure [S] shows the major axis of the sky position error 
ellipse for the same binaries as Fig. [1] except with higher 
harmonics now added to the waveform model. The me- 
dian value of 2a for random-spin, gas- free systems is 21.7 
arcminutes. Comparing to the leading quadrupole wave- 
form value of 34.8 arcminutes, we see that higher har- 
monics do indeed add some additional information not 
contained in precession. The difference between the two 
values is a factor ~ 1.6, consistent with the results of [l6| . 
For partially aligned systems, the shape of the plot shows 
that the higher harmonics have had an important effect; 
both partially aligned histograms look much more like the 
roughly flat random-spin case than the strongly peaked 
histograms shown in Fig. [TJ The median is 28.2 arcmin- 
utes for hot gas and 32.7 arcminutes for cold gas. Both 
results are great improvements on the leading quadrupole 
waveform values (factors ~ 2.2 and 2.8, respectively), 
indicating that the inclusion of higher harmonics has in- 
deed "made up" for the loss of some spin precession. Both 
results are actually better than the leading quadrupole, 
gas-free result of 34.8 arcminutes. In this case, a full 
waveform with a small amount of precession does bet- 
ter than a leading quadrupole waveform with potentially 
significant precession. 

Figure ITOl shows the results for the minor axis 2b, with 
medians 16.1 arcminutes for random spins and 13.7 ar- 
cminutes for both 30° and 10° alignment. These are all 
better than the leading quadrupole, random-spin result 
of 24.6 arcminutes. More interestingly, we see that when 
higher harmonics are included, the "cold gas" and "hot 
gas" errors are smaller than the "no gas" errors. As dis- 
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mi (Mq) 


m2(M0) 


No 
2a 


gas 
2b 


Hot 
2a 


gas 
2b 


Colo 
2a 


gas 
2b 


10^ 


10^ 


21.8 


12.3 


30.3 


14.2 


35.6 


15.9 


3 X 10^ 


10'' 


14.3 


9.67 


19.9 


10.1 


31.0 


13.2 


3 X 10^ 


3 X 10^ 


26.2 


16.2 


30.4 


18.3 


39.1 


25.6 


10'^ 


10" 


14.0 


11.9 


13.7 


9.25 


17.6 


9.69 


10'' 


3 X 10" 


21.7 


16.1 


28.2 


13.7 


32.7 


13.7 


10'' 


10" 


48.1 


32.0 


60.4 


37.1 


53.0 


32.4 


3 X 10^ 


3 X 10" 


29.1 


25.8 


25.7 


20.0 


35.8 


24.4 


3 X 10^ 


10" 


36.0 


26.8 


48.2 


27.4 


58.2 


30.2 


3 X 10" 


3 X 10" 


63.5 


39.8 


103 


58.1 


109 


66.5 


10^ 


10" 


36.7 


32.4 


38.9 


27.4 


54.7 


31.3 


10^ 


3 X lO" 


45.0 


32.8 


65.1 


33.0 


82.7 


36.5 


10^ 


10^ 


114 


65.6 


144 


80.1 


228 


115 



TABLE V: Median sky position major axis 2a and minor axis 
2b, in arcminutes, for binaries of various masses at z — 1, in 
the "no gas" (random-spin), "hot gas" (30° alignment), and 
"cold gas" (10° alignment) cases when higher harmonics are 
included in the waveform model. Bold entries are those that 
do better than the no gas case when higher harmonics are 
ignored (i.e. Table H}. Italic entries do worse than that case, 
but only by 10% or less. 



cussed in the previous section, this is due to the improve- 
ment in SNR in the aligned case. Higher harmonics break 
degeneracies well enough that it is more beneficial to have 
partially aligned binaries with more SNR and less preces- 
sion than randomly aligned binaries with less SNR and 
more precession. 

Table fVl shows results for a variety of masses. All show 
improvement from Table HI however, the improvement 
is not always as strong as in the case discussed above 
(toi = 10^ Mq,TO2 = 3 X 10^ Mq). Bold text indicates 
cases in which the errors match or improve upon the re- 
sults from the leading quadrupole waveform for random 
spins. Italics indicate errors which are worse, but by no 
more than 10%. Because of statistical issues, these cases 
could very well be "bold" in a different Monte Carlo run, 
so we will consider them as such for purposes of summa- 
rizing the results. While hot gas (30° alignment) systems 
achieve this particular benchmark for a majority of mass 
cases, cold gas (10° ahgnment) systems do not. Cold 
gas systems do, however, meet it for a majority of mass 
cases if only the minor axis is considered; 2b generally 
fares better than 2a. Both axes exhibit cases where errors 
are smaller with alignment than without, the minimum 
sometimes occurring in a "sweet spot" of 30° alignment 
and sometimes at 10° alignment. In general, errors are 
better for larger mass ratios. This is to be expected be- 
cause both higher harmonics and precession have a more 
complicated structure for larger mass ratios. Finally, the 
improvements are worst for the smallest masses, where 
the higher harmonics (except the cos$oib terms, which 
are technically "lower" harmonics) begin to go out of 
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FIG. 11: Same as Fig. [9l but for the fractional error in lumi- 
nosity distance, A.Dl/Dl. 




^Qr^ 

Ami/ mi 



10" 



10 



FIG. 12: Same as Fig. [9l but for the fractional error in mass, 
Ami/mi. 



band. 

Figure [11] shows the results for luminosity distance 
errors IS.Dl/Dl. Here the medians are 3.20 x 10^^, 
3.20 X 10'^, and 3.54 x 10~^ for the no gas, hot gas, 
and cold gas cases, respectively. Again, these are all 
better than the leading quadrupole, no gas value of 
5.24 X 10"'^, as we might expect since distance deter- 
mination is strongly tied to sky position determination. 
Table |VT] gives the results for various masses. Most cases 
beat the leading quadrupole, no gas values of Table [Tl] 
or come within 10%. That is, except for the lowest mass 
systems, using the full waveform essentially always brings 
the distance errors for aligned spins back to the level of 
random spins. In this respect, distance errors are similar 
to (and even a bit better than) the minor axis of the sky 
position error ellipse. Finally, as with sky position, some 
mass cases feature errors which decrease as spins become 



aligned. 

Figs. [12] and [13] show the results for masses mi and 
TO2. The medians for Ami/mi are 1.44 x lO^'^ (no gas), 
1.17 X 10"^ (hot gas), and 1.24 x 10"^ (cold gas). For 
Am2/m2, these are 1.15 x 10"^ (no gas), 9.32 x lO"* (hot 
gas), and 9.85 x lO""* (cold gas). These are all significant 
improvements on the leading quadrupole, no gas case, 
about a factor of 3-4. It seems that for mass errors, 
higher harmonics are more useful than spin precession. 

This conclusion is supported by Table [Vlll Here we see 
that every case is significantly improved over the leading 
quadrupole, random-spin result. The effect is strongest 
at higher masses, where the improvement can be an or- 
der of magnitude or more. This is due to the well-known 
effect of higher harmonics on higher mass signals: Nor- 
mally these signals are only in band for a short amount of 
time. The inclusion of higher frequencies keeps the sig- 



mi (Afo) 


ma (Mo) 


No gas 


Hot gas 


Cold gas 




10^ 


3.83 X 10"^ 


5.95 X IQ-^ 


7.23 X 10"^ 


3 X 10^ 


lO'^ 


1.89 X 10"^ 


2.88 X 10"^ 


4.05 X 10"^ 


3 X lO'^ 


3 X 10^ 


4.16 X 10"^ 


4.78 X 10"^ 


5.45 X 10"^ 


10** 


lO'^ 


2.07 X 10"^ 


1.66 X 10"^ 


2.12 X 10"^ 


10'^ 


3 X lO-"' 


3.20 X 10"^ 


3.20 X 10"^ 


3.54 X 10"^ 


10'' 


lO'' 


7.16 X 10"^ 


7.62 X 10"^ 


7.48 X 10~^ 


3 X lO'' 


3 X 10^ 


4.01 X 10"^ 


3.29 X 10"^ 


3.77 X 10"^ 


3 X lO'' 


10" 


5.38 X 10"^ 


6.23 X 10"^ 


6.79 X 10"^ 


3 X 10'^ 


3 X 10" 


0.0115 


0.0139 


0.0152 


10^ 


10" 


5.69 X 10"^ 


5.69 X 10"^ 


7.42 X 10"^ 


10^ 


3 X 10" 


7.69 X 10"^ 


9.02 X 10"^ 


0.0121 


10^ 


10^ 


0.0223 


0.0270 


0.0302 



TABLE VI: Same as Table (V] but for the fractional error in 
luminosity distance, /S.Dl/Dl- 
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FIG. 13: Same as Fig. [9l but for the fractional error in mass, 
Am2/m2. 
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mi (M©) 


m2(M0) 


No 
Ami/mi 


gas 
Am2/m2 


Hot 

Ami/mi 


gas 
Am2/m2 


Cold 
Ami /mi 


gas 
Am2/m2 


10^ 


10'' 


1.22 X 10^ 


1.21 X 10^ 


1.16 X 10"^ 


1.16 X 10"^ 


1.24 X 10"^ 


1.24 X 10"^ 


3 X 10^ 


10^ 


1.29 X 10"^ 


1.04 X 10"^ 


1.43 X 10"^ 


1.16 X 10"^ 


1.77 X 10"^ 


1.44 X 10"^ 


3 X 10^ 


3 X lO'"' 


6.14 X 10 ■* 


6.13 X 10 ■* 


5.37 X lO"'' 


5.35 X 10"" 


5.12 X 10 " 


5.14 X 10 " 


10^ 


10^ 


1.32 X 10^ 


9.19 X 10 ■* 


9.23 X 10 " 


6.39 X 10 


1.10 X 10^ 


7.64 X 10""* 


10" 


3 X lO'"' 


1.44 X 10"^ 


1.15 X 10^ 


1.17 X 10"^ 


9.32 X 10 " 


1.24 X 10"^ 


9.85 X 10"" 


lO'' 


10" 


1.27 X 10"^ 


1.28 X 10"^ 


7.64 X lO"'' 


7.65 X 10"" 


7.25 X lO"" 


7.21 X 10"" 


3 X 10'' 


3 X 10" 


2.50 X 10^ 


1.73 X 10"^ 


1.47 X 10"^ 


1.00 X 10^ 


1.62 X 10"^ 


1.12 X 10"^ 


3 X 10^ 


10" 


3.06 X 10^ 


2.45 X 10"^ 


2.55 X 10"^ 


2.05 X 10"^ 


2.65 X 10^ 


2.13 X 10^ 


3 X 10® 


3 X 10" 


2.04 X 10^ 


2.03 X 10^ 


1.48 X 10"^ 


1.48 X 10"^ 


1.47 X 10"^ 


1.47 X 10"^ 


10^ 


10" 


4.52 X 10"^ 


3.04 X 10^ 


3.17 X 10"^ 


2.08 X 10"^ 


3.55 X 10^ 


2.35 X 10^ 


10^ 


3 X 10" 


4.04 X 10^ 


3.12 X 10^ 


3.56 X 10"^ 


2.73 X 10"^ 


4.21 X 10"^ 


3.28 X 10"^ 


10^ 


10^ 


5.05 X 10^ 


4.90 X 10^ 


4.17 X 10"^ 


4.01 X 10"^ 


4.50 X 10^ 


4.18 X 10"^ 



TABLE VII: Same as Table El but for the mass errors Ami /mi and Am2/m2. 



nal in band longer, allowing for the accumulation of more 
phase and better mass determination. The improvement 
is also greater for equal masses, similar to what was seen 
in Sec. IIIII for precession. 

While the extrinsic parameter errors were only occa- 
sionally reduced by partial alignment, this phenomenon 
occurs almost always for mass errors. This might be ex- 
pected, since the effect showed up previously for mass er- 
rors (in Table nil)) even without higher harmonics to help 
break degeneracies. In general, the difference in mass 
accuracy between gas environments is relatively small. 
Partial alignment of spins does not affect mass determi- 
nation as long as the signal model includes higher har- 
monics; such alignment may even help measure mass, at 
least slightly. 

Finally, Figs. [U and [15] present the errors in spin mag- 



nitude. These figures clearly show that higher harmonics 
do not help spin errors as much as spin precession does. 
This is to be expected, as the spin magnitudes drive the 
precession but do not appear in the higher harmonic am- 
plitudes. Any gain in spin accuracy due to higher har- 
monics is a result of improvement in other parameters 
(such as the masses) which are correlated with the spins. 
For xij the median errors are 2.20 x 10~^ for no gas, 
3.01 x 10"3 for hot gas, and 6.88 x 10"^ for cold gas; 
for X2, these numbers are 1.21 x 10~^, 1.45 x 10~^, and 
3.33 x 10~^. While these errors represent improvements 
(up to a factor of ~ 2) over their leading quadrupole 
waveform counterparts, it is worth noting that for this 
"fiducial" case, spins are the only parameters for which 
the cold gas errors with higher harmonics do not im- 
prove upon the no gas errors for the leading quadrupole 
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FIG. 14: Same as Fig. [9] but for the error in spin magnitude, 
Xi = |Si|/m?. 



FIG. 15: Same as Fig. [9] but for the error in spin magnitude, 
X2 = |S2|/mi. 
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No 

Axi 


gas 

Ax2 


Hot ga 
Axi 


s 

Ax2 


Cold g 
Axi 


IS 

AX2 




lO'^ 


0.0155 


0.0157 


0.0178 


0.0178 


0.0286 


0.0280 


3 X lO'^ 


10^ 


1.79 X 10"^ 


6.18 X 10^ 


3.09 X 10^ 


0.0118 


7.03 X 10"^ 


0.0275 


3 X lO'^ 


3 X 10^ 


0.0207 


0.0211 


0.0247 


0.0243 


0.0335 


0.0335 


10^ 


10^ 


6.43 X 10 " 


0.0196 


6.21 X 10 * 


0.0137 


1.61 X 10"^ 


0.0263 


10'' 


3 X 10^ 


2.20 X 10"^ 


0.0121 


3.01 X 10^ 


0.0145 


6.88 X 10~^ 


0.0333 


10** 


lO'' 


0.0311 


0.0326 


0.0385 


0.0396 


0.0603 


0.0590 


3 X lO'' 


3 X lO'^ 


1.07 X 10"^ 


0.0373 


8.64 X 10 * 


0.0218 


2.24 X 10"^ 


0.0408 


3 X lO'' 


10** 


3.90 X 10^ 


0.0153 


5.29 X 10"^ 


0.0221 


0.0129 


0.0551 


3 X lO'' 


3 X lO'' 


0.0664 


0.0669 


0.0932 


0.0899 


0.172 


0.171 


10^ 


lO'' 


1.75 X 10"^ 


0.0552 


1.67 X 10"^ 


0.0482 


4.03 X 10^ 


0.0872 


10^ 


3 X lO*' 


6.03 X 10^ 


0.0295 


0.0111 


0.0557 


0.0264 


0.140 


10^ 


10^ 


0.495 


0.525 


0.548 


0.581 


1.08 


1.08 



TABLE VIII: Same as Table El but for the spin magnitude errors Axi and Ax2. 
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FIG. 16: Factors by which measurement accuracy improves for different parameters when various degeneracy-breaking effects 
are included in the signal: spin precession (SP), higher harmonics (HH), and both (SP+HH). We also show the product of 
the individual precession and harmonic improvements (SP*HH); this represents the naive limit by which the two effects would 
improve measurement accuracy if their individual improvements simply combined. Each point represents one of the 12 mass 
cases, arranged in order from left to right as they read top to bottom in Tables [II iHl IIIIl llVl IVl IVll IVIII and lVIlH 
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waveform. (Recall, however, that for other masses, sky 
position and distance errors do not always achieve this 
benchmark.) Table IVIIII shows the results for various 
masses. The cold gas error only beats the no gas, leading 
quadrupole error in a few cases. 

We conclude this section by looking more directly at 
the impact of different degeneracy-breaking effects on pa- 
rameter errors. Figure [TH] shows multiplicative improve- 
ment factors (i.e., ratios of errors) when either precession, 
harmonics, or both effects are included in the waveform. 
For this purpose, we consider binaries with 10° alignment 
to represent "no precession" and those with random spins 
to represent "precession." (Remember, though, that even 
10° alignment permits enough precession to significantly 
impact mass and spin estimation.) We also show the 
product of the individual improvements from precession 
and harmonics. This naive limit describes how the im- 
provements would combine if each effect were completely 
uncorrelated from the other. 

For extrinsic parameters (represented in the figure by 
2a; results for 2b and are very similar), the re- 
sults confirm what was known previously. While includ- 
ing both precession and harmonics improves errors more 
than one effect alone, the total improvement falls well 
short of the naive expectation. In essence, a degener- 
acy can only be broken once. For mass (shown on a 
log scale), the same is true, but with the special feature 
pointed out above: Once harmonics are included, they 
essentially dominate mass accuracy determination. Spin 
precession is then a small liability, because the associ- 
ated misaligned spins reduce SNR. For spin errors, on 
the other hand, the combined improvement due to both 
effects does roughly match the naive expectation. In the 
case of xi i the improvement is actually greater than the 
product of the individual improvements. Different be- 
havior for spin magnitudes is to be expected, since the 
information about spin contained in harmonics is very in- 
direct. This information is therefore independent of any 
information derived directly from precession. 



V. CONCLUSIONS 

Past work has shown that both spin precession and 
higher harmonics improve LISA's ability to measure 
the parameters of merging massive black hole binaries. 
Though these two effects produce similar degeneracy- 
breaking effects and similar improvements to measure- 
ment errors, there is one key difference between the two: 
Higher harmonics are always present in the signal (al- 
though the strength of odd harmonics depends on mass 
ratio) . Spin precession, on the other hand, may be highly 
attenuated for physical reasons, namely, the partial align- 
ment of spins due to interaction with gas. In this paper, 
we have studied how this partial alignment affects pa- 
rameter measurement errors. 

Initially ignoring the impact of higher harmonics, we 
found that sky position and distance are measured a fac- 



tor of ~ 1.5 — 2 less accurately for systems aligned within 
30° (due to hot gas) and ~ 2 — 3 less accurately for sys- 
tems aligned within 10° (due to cold gas). A degradation 
of ~ 3 would correspond to an order of magnitude de- 
crease in LISA'S ability to localize a source on the sky and 
a half order of magnitude decrease in the ability to local- 
ize it in redshift space. Since systems with gas, whether 
hot or cold, are the most likely to produce electromag- 
netic counterparts, this means that the results of Paper I 
and [H,[lH strongly overestimate our ability to find these 
counterparts. Mass and spin measurements are also de- 
graded by spin alignment, in some cases by factors up to 
^ 9. However, because the masses and spins are already 
measured quite well, the degradation is not as harmful. 

Adding higher harmonics to the signal substantially 
improves these results. In some cases, measurement er- 
rors for aligned systems can be brought below the error 
level for random spin orientations without higher har- 
monics. For the mass measurements, this improvement 
happens in every case. The minor axis of the sky position 
error ellipse and luminosity distance achieve this bench- 
mark less often, though still in the majority of cases, 
while the major axis and the spin magnitudes do not fare 
as well. Sometimes parameters are actually determined a 
bit better for aligned spins than for random ones thanks 
to the increased SNR measured for aligned spin systems. 
For mass measurements, this happens in almost every 
case we considered. 

Although studies like [l^ and this one are starting 
to finalize expectations for LISA's parameter estimation 
capabilities, several avenues of research into the problem 
still remain. First is the issue of including proper astro- 
physical information, such as the possibility of partially 
aligned spins, when analyzing LISA science. Another ex- 
ample is the recent realization that eccentricity may need 
to be included in the waveform model (with a first anal- 
ysis by Key and Cornish [l^). In order to make reliable 
estimates of LISA's science capabilities, future studies 
must continue to incorporate the newest astrophysical 
developments. It will also be useful to turn the problem 
around and ask what LISA measurements of properties 
like spin alignment can tell us about the surrounding gas 
(or lack thereof). It may be possible, from the gravi- 
tational waves alone, to predict the nature of an elec- 
tromagnetic counterpart or to make statements about a 
binary's environment in case a counterpart is missed. We 
plan to study this issue in more detail in the future. 

Lacking from our Fisher-matrix-based study is an in- 
vestigation of how the spin- alignment priors of the hot, 
cold, and dry scenarios may affect parameter estima- 
tion. In a more thorough Bayesian analysis, the range of 
spin alignments for each model can be included as priors. 
Bayesian model selection can then be used to identify the 
model that best describes the data. For cold gas merg- 
ers, the tight priors on the spin-orbit alignment would 
translate into improvements in the parameter estimation 
over what we have found here. Thus, the spin-alignment 
priors can help put back some of parameter recovery ac- 
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curacy that is taken away by the suppression of the spin 
precession. We are currently investigating this issue. 

An obvious avenue for LISA parameter estimation 
studies is the improvement of the waveform model. While 
the inspiral is nearly complete, the study of the merger 
is just beginning. In the next few years, complete inspi- 
ral waveforms from codes like this one will be joined to 
effective-one-body waveforms for the late inspiral, fits to 
numerical models of the merger, and perturbative ring- 
downs to give a complete LISA waveform for analysis 
purposes. While such studies have begun [13, [l3> ^hey 
do not yet include spins. Including the merger and ring- 
down is critical; not only do they provide a great deal of 
SNR and parameter information, but they also provide 
a physical tapering of the waveform. The results of this 
paper unfortunately do not always quantitatively match 
those of Paper I and [IQ] when appropriate, primarily 
because the earlier studies use the stationary phase ap- 
proximation while we taper the signal and apply an FFT. 
In essence, despite the use of the MECO, we are applying 
an earlier cutoff and losing some information about our 
parameters. When the complete signal is used, choices of 
cutoff will become irrelevant, and different results should 
agree more readily. Of course, for this paper, exact error 
estimates are not the end goal; instead, we have aimed 
only to show general behavior and trends which are in- 
dependent of any shift in the baseline error values. 

Finally, the Fisher-matrix formalism itself must be 
checked to make sure it is correct with complicated wave- 
forms and large numbers of parameters. We plan to 



carry out a comparison between this code's results and 
those obtained by exploring the full posterior probabil- 
ity using Markov Chain Monte Carlo techniques [48l - [5l| . 
Early results have shown that the Fisher matrix is indeed 
still valid in most regimes — a relief given that Markov 
Chain Monte Carlo techniques cannot as easily survey 
a plethora of sky locations and orientations — but the 
full parameter space has not been explored. It is clear, 
however, that care must be taken with all three areas, 
astrophysics, general relativity, and statistical analysis, 
before a final picture of LISA science capabilities can be 
established. 
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